Quantum Griffiths effects in itinerant Heisenberg magnets 



Thomas Vojta 

Department of Physics, University of Missouri-Rolla, Rolla, MO 65409 

Jorg SchmaHan 

Department of Physics and Astronomy and Ames Laboratory, Iowa State University, Ames, I A 50011 

(Dated: February 2, 2008) 

We study the influence of quenched disorder on quantum phase transitions in itinerant magnets 
with Heisenberg spin symmetry, paying particular attention to rare disorder fluctuations. In contrast 
to the Ising case where the overdamping suppresses the tunneling of the rare regions, the Heisenberg 
system displays strong power-law quantum Griffiths singularities in the vicinity of the quantum 
critical point. We discuss these phenomena based on general scaling arguments, and we illustrate 
them by an explicit calculation for 0{N) spin symmetry in the large- Af limit. We also discuss broad 
implications for the classification of quantum phase transitions in the presence of quenched disorder. 



I. INTRODUCTION 

At low temperatures, strongly correlated materials can 
display a surprising sensitivity to small amounts of im- 
perfections and disorder. This effect is particularly pro- 
nounced close to a quantum phase transition, where 
large fluctuations in space and time become fundamen- 
tally connected. Thermodynamic and response prop- 
erties of a material are then affected much more dra- 
matically than close to a classical phase transition, per- 
mitting exotic phenomena like infinite-randomness criti- 
cal points with activated rather than power-law dynam- 
ical scalingSiiiiSiSiLSiS, smeared transitions^^, or non- 
universal exponents at certain impurity quantum phase 
transitions^. 

One interesting aspect of phase transitions in disor- 
dered systems are the Griffiths effectsi^. They are caused 
by large spatial regions that are devoid of impurities and 
can show local order even if the bulk system is in the dis- 
ordered phase. The fluctuations of these regions are very 
slow because they require changing the order parameter 
in a large volume. GrifBthsi^ showed that this leads to 
a singular free energy in a whole parameter region in the 
vicinity of the critical point which is now known as the 
Griffiths phase. In generic classical systems, the contri- 
bution of the rare regions to thermodynamic observables 
is very weak since the singularity in the free energy is 
only an essential oneiSiiii^. The consequences for the 
dynamics are more severe with the rare regions dominat- 
ing the long-time behaviorii^iiSii&. 

Due to the perfect disorder correlations in (imagi- 
nary) time, Griffiths phenomena at zero temperature 
quantum phase transitions are enhanced compared to 
their classical counterparts. In random quantum Ising 
systemsSiSi^i^ and quantum Ising spin glasse o'^i^i^ , ther- 
modynamic quantities display power-law singularities 
with continuously varying exponents in the Griffiths 
phase, with the average susceptibility actually diverging 
inside this region. 

The systems in which these quantum Griffiths phe- 
nomena have been shown unambiguously all have un- 



damped dynamics (a dynamical exponent z = 1 in the 
corresponding clean system). However, many systems 
of experimental importance involve magnetioiiiSiiSiSS or 
superconducting^^ degrees of freedom coupled to conduc- 
tion electrons which leads to overdamped dynamics char- 
acterized by a clean dynamical exponent z > 1. Study- 
ing the effects of rare regions in this case is therefore an 
important issue. In recent years, there has been an in- 
tense debate on the theory of quantum Griffiths effects 
in itinerant Ising magnets. It has been suggested^ that 
overdamped systems show quantum Griffiths phenomena 
similar to that of undamped systems. However, recently 
it has been showniS*2^ that the overdamping prevents the 
rare regions from tunneling leading to static rare regions 
displaying superparamagnetic rather than quantum Grif- 
fiths behavior, at least for sufficiently low temperatures. 
In Refi-^ it was also pointed out that different behavior 
is expected for systems with continuous spin symmetry. 

In this paper, we examine quantum Griffiths effects 
in itinerant Heisenberg magnets. Our results can be 
summarized as follows: In contrast to the Ising case, 
itinerant magnets with Heisenberg symmetry (in gen- 
eral, O (N) symmetry with iV > I) do display power-law 
quantum Griffiths singularities. The locally ordered rare 
regions are not static but retain their quantum dynam- 
ics. Their low-energy density of states follows a power 
law, /9(e) ^ e''/^ where d is the space dimensional- 
ity and z' is a continuously varying dynamical exponent. 
This leads to power-law dependencies of several observ- 
ables on the temperature T, including the specific heat, 
C ~ T'^/^ , and the magnetic susceptibihty, x T^^^ 
Our results are not limited to Heisenberg magnets, they 
generally apply to O {N) order parameters with > I 
including the superconductor-metal transition^^ in thin 
nano-wireaSi. 

The paper is organized as follows: In section ^ we 
derive quantum Griffiths effects from general scaling ar- 
guments based on the observation that a rare region in 
an itinerant Heisenberg magnet is at its lower critical 
dimension. These arguments suggest a general classifica- 
tion of disordered quantum phase transitions in terms of 
the dimensionality of the rare regions. In section UTTl we 
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then present an explicit calculation for 0{N) spin sym- 
metry in the large- limit. We conclude in section lTVI bv 
discussing the importance of the spin symmetry, the re- 
lation to Kondo physics as well as possible experimental 
realizations of our predictions. 



II. QUANTUM GRIFFITHS EFFECTS FROM 
SCALING ARGUMENTS 

Our starting point is a quantum Landau-Ginzburg- 
Wilson free energy functional for an A^-component {N > 
1) order parameter field (j) = ■ • • j 07v)- For definite- 
ness, we consider the itinerant antiferromagnetic transi- 
tion. The action of the clean system reads2S*2Si2i 

Jdxdy 0(x) Fix, v) m + ^Jdx (t>\x) . (1) 

Here, x = (x, r) comprises position x and imaginary time 

r, and J dx = / dx /q^''^ dr. The imaginary time direc- 
tion which characterizes the quantum dynamics formally 
appears like an additional spatial dimension of a classical 
system. The perfect disorder correlations in this imag- 
inary time direction are ultimately responsible for the 
enhancement of the Griffiths effects at zero temperature. 
T{x,y) is the bare inverse propagator (bare two-point 
vertex), whose Fourier transform is 

r(q,c^„) = (^o + q'+7kn|'/') (2) 

and ro is the bare energy gap, i.e., the bare distance from 
the clean critical point. 

We are interested in the case of overdamped spin dy- 
namics (z = 2) with 7 ~ (Jpp)^ I [EfOq) where J is the 
coupling constant between spin degrees of freedom and 
conduction electrons, with DOS, pp, and Fermi energy, 
Ep, respectively, ag is the lattice constant. In order to 
demonstrate the special behavior of Heisenberg systems 
with z = 2, where overdamping is caused by the particle- 
hole continuum of itinerant electrons, we frequently dis- 
cuss variable z and compare the behavior with ballistic 
spin systems with z = 1. We will use a system of units 
with 7 = 1. The clean system undergoes the quantum 
phase transition when the renormalized gap r vanishes. 
To introduce quenched disorder, we dilute the system 
with nonmagnetic impurities of spatial density p, i.e., we 
add a random potential, (5r(x) — J2i ~ to tq. 

Here, x(i) are the the positions of the impurities, and 
V^(x) is a positive short-ranged impurity potential. 

We first present the general scaling arguments leading 
to quantum Griffiths behavior in this system. Despite the 
dilution, there are statistically rare large spatial regions 
devoid of impurities and thus unaffected by the disorder. 
The probability for finding such a region of volume L'^, 
frequently referred to as an instanton, is 

w - (1 - = exp(-cL'^) (3) 



with c = —a^ ln(l — p). Below the clean critical point, 
the rare regions can be locally in the ordered phase even 
though the bulk system is not. At zero temperature, each 
rare region is equivalent to a one-dimensional classical 
0{N) model in a rod-like geometry: finite in the d space 
dimensions but infinite in imaginary time. For over- 
damped dynamics, z = 2, the interaction in imaginary 
time direction is of the form (t — t')~^. One-dimensional 
continuous-symmetry 0{N) models with 1/x^ interac- 
tion are known to be exactly at their lower critical 
dimensionSSiSSi^. Therefore, an isolated rare region of 
linear size L cannot independently undergo a phase tran- 
sition. Its energy gap depends exponentially on its vol- 
ume (i.e., the effective spin of the droplet) 

eL ~ expi-bL") . (4) 

Equivalently, the susceptibility of such a region diverges 
exponentially with its volume. Combining Q and 
gives a power-law density of states for the energy gap e 
(to leading exponential accuracy), 

p{e) oc e=/''-i = e'^/"'-^ (5) 

where the second equality defines the customarily used 
dynamical exponent z'>^. It continuously varies with dis- 
order strength or distance from the clean critical point. 
Many results follow from this. For instance, a region with 
a local energy gap e has a local spin susceptibility that 
decays exponentially in imaginary time, xioc(''' ~> oo) cx 
exp(— er). Averaging by means of p yields 

Xr:c(^^oo)(xr-'^/^'. (6) 

The temperature dependence of the static average sus- 
ceptibility is then 

Xioc(T) = dr x^cir) « T'/''-\ (7) 

Jo 

If d < z' , the local zero-temperature susceptibility di- 
verges, even though the system is globally still in the 
disordered phase. Analogously, the contribution of the 
rare regions to the specific heat C can be obtained from 

AE = J de pie) e e~^/^/(l -t- e"^/^) cx T'^Z-'+i (8) 

which gives AC cx T'^^^ . Other observables can be de- 
termined in a similar fashion. The power-law density of 
states in the Griffiths phase of a disordered itinerant 
0{N) magnet and the resulting quantum Griffiths singu- 
larities lO, 0, l|Hll are the central results of this Paper. 
They take the same form as the quantum Griffiths sin- 
gularities in undamped (clean z = 1) random quantum 
Ising modelsSi^i^i^ and quantum Ising spin glassesi^. 

These scaling arguments suggest a general classifi- 
cation of Griffiths phenomena in the vicinity of bulk 
phase transitions (at least those described by Landau- 
Ginzburg- Wilson theories) with weak, random-T^ or 
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random-mass type disorder. It is based on the effective 
dimensionality of the rare regions. Three cases can be 
distinguished: 

(i) If the rare regions are below the lower critical di- 
mensionality of the problem, their energy gap depends 
on their size via a power law, ^ L""^. Since the prob- 
ability for finding a rare region is exponentially small in 
L, the low-energy density of states in this first case is 
exponentially small. This leads to weak "classical" Grif- 
fiths singularities characterized by an essential singular- 
ity in the free energy. This case is realized in generic 
classical systems (where the rare regions are finite in all 
directions and thus effectively zero-dimensional). It also 
occurs in quantum rotor systems with Heisenberg sym- 
metry and undamped (z — 1) dynamice'^^'^^. Here, the 
rare regions are equivalent to one-dimensional classical 
Heisenberg models which are also below = 2. 

(ii) In the second class, the rare regions are exactly 
at the lower critical dimension. In this case, their en- 
ergy gap shows an exponential dependence [like (@}] on 
L. As shown above, this leads to a power-law density 
of states and strong power-law quantum Griffiths singu- 
larities. This second case is realized, e.g., in classical 
Ising models with linear defects2i^ and random quan- 
tum Ising models (each rare regions corresponds to a 
one-dimensional classical Ising model)^' as well as in the 
disordered itinerant quantum Heisenberg magnets stud- 
ied here (the rare regions are equivalent to classical one- 
dimensional Heisenberg models with interaction). 

(iii) Finally, in the third class, the rare regions are 
above the lower critical dimension, i.e., they can undergo 
the phase transition independently from the bulk sys- 
tem. In this case, the locally ordered rare regions be- 
come truly static which leads to a smeared phase tran- 
sition. This happens, e.g., for classical Ising magnets 
with planar defects^ (the rare regions are effectively two- 
dimensional) or for itinerant quantum Ising magnetsiSiS^ 
where the rare regions are equivalent to classical one- 
dimensional Ising models with interaction. 



III. RARE REGIONS IN THE LARGE-iV LIMIT 

To complement the general scaling arguments and to 
obtain quantitative estimates for the exponent z' we now 
perform an explicit calculation of the Griffiths effects in 
the model ^ in the large- limit. The approach is a 
generalization of Bray's treatmen'ri^ of the classical case. 
In the large- A^ limit, a clean system undergoes a quan- 
tum phase transition for g = gc (x A^~<^~^ with coupling 
constant g — and upper momentum cut off A. For 
g < gc, the clean system is in the ordered state with the 
order parameter 



rO, clean 



[^(5c-<?)/(<?c5)]^/^ 



(9) 



and vanishing gap. In the random system, we consider 
a droplet of size L'^, devoid of impurities and determine 



its size dependent energy gap, e. It is determined by the 
equation of state e(f>o = h, where 



N ' 



(10) 



h is the field conjugate to the order parameter, (j)o = (0) , 
of the droplet and 



TL- 



,2/z- 



(11) 



For T > 0, both the q and u! sums are discrete. Con- 
sequently, the order parameter 0o — h/e vanishes for 
h ^ since e must remain positive to avoid a diver- 
gence of the q = 0, w„ = contribution to (0^). Thus, 
classical droplets are below d~ . t T = 0, a frequency 
integration must be performed and the e ^ limit be- 
comes less singular. Yet, for z < 2 droplets remain below 
d~ since the q = contribution to (0^) still diverges as 

L~'^e^~ . For z = 2 this term diverges only as In (eL^) 
and, as expected, droplets with z = 2 are marginal and 
located at their lower critical dimension. 

To quantify these arguments and to determine the de- 
pendence of e on L for T = 0, we apply the finite size 
analysis of the large- A^ theory"^^ to the quantum limit. 
As shown in the appendix, we obtain for eL^ <C 1 and 
z = 2: 



1 L 



-d 



h^) = ±-— In (eL^) . 

5c TT 



(12) 



Inserting this into H1U|) gives for small e: 

e = L-^ c^p {-bL"^) , (13) 



with b 



9cS 



N 



'o clean- This explicitly verifies eq. 
in the large-A^ limit. For a given distance, b, to the 
clean critical point, only droplets larger than a certain 
value contribute to quantum Griffiths behavior. For z < 
2, the last term in (|12|l is proportional to L^'^e^~ and 
we obtain e cx with -0 — ^r^- For z > 2, (j)o ^ 0, i.e. 
even a finite droplet is allowed to order at T = 0. Since 
the onset of order depends on the size of the droplet a 
smearing of the transition occurs. All this is in agreement 
with our general expectation, discussed above. 

Inserting our result for the coefficient b into we 
obtain an explicit expression for the Griffiths exponent 



rfTT {gc - g) Qq 
5cgln(l -_p)-i 



(14) 



The density of states (jSJ diverges for e ^ if z' > d. 
z' vanishes as one approaches the clean critical point 
g ^ gc, but becomes larger as {gc — g) /g grows. In Fig. 
1, we plot the coupling constant g* , below which z' > d, 
as function of the impurity concentration, p, for three 
different values of the non-universal number agA. Quan- 
tum Griffiths effects dominate the low energy excitations 
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FIG. 1: Coupling constant g* / Qc below which quantum Grif- 
fiths effects cause a diverging low energy density of states, as 
function of disorder concentration, p, in two and three dimen- 



sions for various values of A = 



for g < g*, provided that droplets are still sufhciently 
diluted and the system is above the critical point, g^^^, of 
the random system. Observable quantum Grifhths effects 
exist unless Aoq becomes small. 

At finite temperatures, a crossover occurs to weaker 
classical Griffiths effects. To estimate the characteristic 
crossover temperature for z = 2, we decompose ((^^), eq. 

into its zero-temperature part and the more singular 
classical (a;„ = q = 0) contribution: 



-d 



T 
El ' 



(15) 



The crossover occurs when the classical term becomes 
comparable to {^(jy^'j We find that droplets with L > 

Lo (T), determined by T = (6/7r)Lg~^e"''^o , behave clas- 
sically and p (e) is suppressed for e < = LQ^e~^^o . 
Droplets smaller than Lq (T) still follow the quantum dy- 
namics. Quantum Grifhths behavior persists as long as 
Co (T) < T. This is fulfilled for sufficiently low tempera- 
tures T < To = fdb'^''^ with fd = 7r-2/'iexp(-7r). Above 
the temperature Tq the density of states for w > T is 
then given as 



Pciass (e) oc exp 



A 



(16) 



with A = ^ ln(l — p) ^, in agreement with the behavior 

for classical Heisenberg systems^-. Instead of adding the 
purely classical contribution to (^(fP')^_^, the behavior at 
low T may be described by explicitly calculating the first 
low temperature corrections to (^(fp'^^_^. We find that 
these low temperature corrections behave as 



L- 



J.2 



487r {4neL) 



(17) 



An argumentation identical to the one given below Ea ll5l 
yields that these corrections become relevant above a 



temperature which behaves identical to Tq cx 6^/'', only 
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with a different numerical prefactor fd — (l67r\/3) . 
The closer one one approaches the clean quantum phase 
transition, where b vanishes, the narrower is the region of 
quantum Griffiths behavior. More importanly, very close 
to the dirty critical point, g^^^, the crossover temperature 
is exponentially suppressed because the droplets have a 
minimum size of the order of the bulk correlation length. 

IV. DISCUSSION AND CONCLUSIONS 

To summarize, we have studied quantum Griffiths ef- 
fects in itinerant magnets with continuous order parame- 
ter symmetry, using the itinerant antiferromagnet as the 
primary example. We have shown that this system dis- 
plays strong power-law quantum Griffiths singularities. 
In this section we will address a few open questions and 
important implications of the results. 

Our first point concerns the importance of Kondo 
physics. In eq. 1^ we assumed, following RefsiS^iSSiSi, 
that spin degrees of freedom in disordered metals can 
be described by quantum rotors with overdamped dy- 
namics. One might worry about the Kondo dynamics of 
the entire dropletwhich is not included in the rotor ap- 
proach (in Ref^ it was shown that an extended mag- 
netic structure in a metallic environment can behave 
at low T as an ansiotropic multi-channel Kondo prob- 
lem). Our theory is valid only if the fc-channel Kondo 
behavior, presumably with large k ~ (L/ao) , emerges 
only for Tk <C £ (L) . Using the results for the related 
problem of a large droplet induced by a single magnetic 
impuritySSi^ it indeed follows that the Kondo tempera- 
ture Tfc ~ e (i) exp [— const. L'^^ is exponentially smaller 
than the crossover scale for quantum Griffiths behavior 
and thus negligible. 

We emphasize the difference between continuous spin 
symmetry and Ising symmetry. For Ising symmetry, rare 
regions are above the lower critical dimension. They 
cease to tunnel and become static at sufficiently low tem- 
peratures, leading to superparamagnetic behavior^^ and, 
ultimately, to a smeared transitionlS. Quantum Griffiths 
behavior can at best occur in a transient temperature 
window^S. In contrast, for continuous symmetry, the 
rare regions are exactly at the lower critical dimension 
and retain their dynamics, with a power- law low-energy 
density of states. Quantum Griffiths effects dominate the 
low-temperature physics for g* > g > g^'^. 

Griffiths phenomena in itinerant /erromagnets require 
separate attention because mode-coupling effects induce 
a long-range interaction of the spin fluctuations^". This 
can potentially change the conditions for locally ordered 
droplets and thus the form of the Grifhths effects. 

Let us comment on measuring the predicted effects. 
Many of the heavy electron systems displaying magnetic 
quantum phase transitions have a strong spin anisotropy 
and are thus better described by Ising models. Gd-based 
intermetallics have a local Heisenberg symmetry, but the 
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hybridization between magnetic and conduction electrons as well as 
is very small. This yields a very low temperature, Tg, for 
the onset of quantum Griffiths behavior. Most promis- 
ing are Sd-transition metal systems with weak spin-orbit 
interaction and strong hybridization. A candidate is the with 
3d heavy fermion system LiV20^i where recent experi- 
ments did show a broad distribution of relaxation rates^. 
In addition, the XF-version of our theory {N = 2) di- 
rectly predicts quantum Griffiths behavior in disordered 
thin nanowiresSi close to the metal-superconductor quan- 
tum phase transitions^, with direct impact on the con- 
ductance fluctuations of these systems. Our results hope- 
fully motivate further the search for new and unconven- 
tional behavior in dissipative quantum systems with con- 
tinuous order parameter symmetry. 

This letter has focused on the Griffiths region above 
the dirty critical point. There is, however, a possible 
connection to the properties of the quantum critical point 
itself. It is known that the quantum critical points of un- 
damped random quantum Ising models, which also dis- 
play power-law quantum Griffiths effects, are of exotic 
infinite-randomness typeSiSiii^. The underlying strong- 
disorder renormalization group^'"*^ supports a close con- 
nection between the quantum Griffiths effects and the ex- 
otic critical properties. This suggests that the quantum 
critical point of disordered itinerant Heisenberg magnets 
may also be of infinite-randomness type. 
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Energy by Iowa State University under Contract No. W- with 
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J=^5(n) 



d'^q f duj e^i "^ 



2/z 



(20) 



(21) 



Here n — (ni, ...rid) is a d-component vector integer com- 
ponents. 

As long as c?-f z > 2, which we assume throughout this 
paper, it holds for small e that {4>'^) ^ — '^here gc is 
the critical coupling constant of the clean bulk quantum 
phase transition, i.e. for L oo. The corrections behave 
as e 2 for d -t- z < 4, as eA'^~^^~^ for d + z > i, and 
as elog (A^/e) for d + z = 4, i.e. vanish as e ^ 0. 

Next we analyze the sum over n in Eql^Ul which takes 
into account the finite size of the droplet. Using ^ — 
daexp (—ax), it follows 



S{n) 



dae-"" 



J (,27r) 



n 



^—\,..,d 



where we performed the momentum and frequency inte- 
grations. Substituting t — 47raL~^ it we obtain for J: 



r (i + L'^-{d+^) 



A{p) = 



V. FINITE SIZE ANALYSIS OF THE LARGE iV 

THEORY and 



dt _£t_ 
e 4" 



t{d+z)/2 



B{-\ -1 



(22) 



(23) 



In this appendix we summarize the derivation of Ea ll2l 
by generalizing the finite size analysis of Ref,"^^-. to the 
quantum case. The order parameter fluctuations 



TL 



-d 



|2/z 



(18) 



determine the value of e {L) through e — ri^+u (0^) + "^- 
We consider a droplet of size L"* and made the assump- 
tion that there are no impurities inside of it. Then we 
are allowed to decsribe this droplet as a clean, but finite 
system. Considering the zero temperature limit, the fre- 
quency summation becomes an integration. The discrete 
sum over momenta is then analyzed using the Poisson 
summation formula and we obtain l^cjp) = {'P'^) ^ + J 
with 



B{t)^ J2 exp(- 



-irtn 



It holds B (t) 



{j)^^^B{\) which leads to 



-Pi 

-e 



dt lB{t) 



- 1 



e 4" 



(24) 



(25) 



For z < 2, the dominant contribution to A (p) comes 
from the first integral in Eg 1251 which can be evaluated 
explicitly. The second integral is finite as p = eL^ 
for all z. For small p and z < 2 follows for the leading 
contributions 



duj 



1 



27r q2 



,2/2 



(19) 



Hp) 



r(i 



(47r) 



- + O {p') 



(26) 
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while z = 2: 

A(p)^-log(£)+0(p") (27) (^^).l-^log(eX^) (28) 

9c TT 

if z = 2. If z > 2, ^ (p 0) remains finite. 

Collecting the various contributions to (0^) yields for which is the result Eall2l given above. 



^ B.M. McCoy, Phys. Rev. Lett. 23, 383 (1969). 
^ D.S. Fisher, Phys. Rev. Lett. 69, 534 (1992); Phys. Rev. 
B 51, 6411 (1995). 

A.P. Young and H. Rieger, Phys. Rev. B 53, 8486 (1996). 
^ C. Pich, A.P. Young, H. Rieger, and N. Kawashima, Phys. 

Rev. Lett. 81, 5916 (1998). 
® O. Motrunich, S.-C. Man, D.A. Huse, and D.S. Fisher, 

Phys. Rev. B 61, 1160 (2000). 

M. ThiU and D. Huse, Physica A 214, 321 (1995). 
* M. Guo, R. Bhatt and D. Huse, Phys. Rev. B 54, 3336 

(1996). 

^ H. Rieger and A.P. Young, Phys. Rev. B 54, 3328 (1996). 
1° T. Vojta, Phys. Rev. Lett. 90, 107202 (2003). 
" A. Georges and A.M. Sengupta, Phys. Rev. Lett. 74, 2808 

(1995), and references therein. 
12 R.B. Griffiths, Phys. Rev. Lett. 23, 17 (1969). 
1^ A.J. Bray and D. Huifang, Phys. Rev. B 40, 6980 (1989). 
^'^ M. Randeria, J. Sethna, and R.G. Palmer Phys. Rev. Lett. 

54, 1321 (1985). 
1^ A.J. Bray, Phys. Rev. Lett. 59, 586 (1987). 
1® D. Dhar, M. Randeria, and J. P. Sethna, Europhys. Lett. 

5, 485 (1988). 

1^ H. V. Lohneysen et al., Phys. Rev. Lett. 72, 3262 (1994). 

1® S.A Grigera et al.. Science 294 329 (2001). 

1^ C. Pfleiderer, G.J. McMuUan, S.R. Juhan, and G.G. Lon- 

zarich, Phys. Rev. B 55, 8330 (1997). 
2° M. C. de Andrade et al., Phys. Rev. Lett. 81, 5620 (1998). 

21 A. Rogachev and A. Bezryadin, Appl. Phys. Lett. 83, 512 
(2003). 

22 A.H. Castro Neto, G. Castilla, and B.A. Jones, Phys. Rev. 
Lett. 81, 3531 (1998); A.H. Castro Neto and B.A. Jones, 
Phys. Rev. B 62, 14975 (2000). 

23 A.J. Millis, D.K. Morr, and J. Schmalian, Phys. Rev. Lett. 



87, 167202 (2001); Phys. Rev. B 66, 174433 (2002). 
2* S. Sachdev, P. Werner, and M. Troyer, Phys. Rev. Lett. 

92, 237003 (2004). 
25 J. Hertz, Phys. Rev. B 14, 1165 (1976). 
2« A.J. Millis, Phys. Rev. B 48, 7183 (1993). 
2^ T.R. Kirkpatrick and D. Belitz, Phys. Rev. Lett. 76, 2571 

(1996); 78, 1197 (1997). 
2* G.S. Joyce, J. Phys. C 2, 1531 (1969). 
2^ F.J. Dyson, Commun. Math. Phys. 12, 91 (1969). 
3° P. Bruno, Phys. Rev. Lett. 87, 137203 (2001). 

31 A.P. Young, Phys. Rev. B 56, 11691 (1997). 

32 N. Read, S. Sachdev and J. Ye, Phys. Rev. B 52, 384 

(1995) . 

33 R. Sknepnek, T. Vojta, and M. Vojta, 'cond-mat/0402352" 
3'' B.M. McCoy and T.T. Wu, Phys. Rev. 176, 631 (1968); 

188, 982 (1969). 
35 T. Vojta, J. Phys. A 36, 10921 (2003); R. Sknepnek and 

T. Vojta, Phys. Rev. B 69, 174410 (2004). 
3'5 E. Brezin, J. Phys. (Paris) 43, 15 (1982). 
3^ N. Shah and A. J. Millis, Phys. Rev. Lett. 91, 147204 

(2003). 

38 A. L Larkin and V. L Melnikov, Sov. Phys. JETP 34, 656 
(1972). 

3^ Y.L. Loh, V. Tripathi, and M. Turlakov, 

cond-mat/0405618 
■1° T.R. Kirkpatrick and D. Belitz, Phys. Rev. B 53, 14364 

(1996) . 

''1 S. Kondo et al, Phys. Rev. Lett. 78, 3729 (1997). 

''2 H. Kaps et al., J. Phys.: Condens. Matter 13, 8497, (2001); 

D.C. Johnston, et aZ. (unpublished). 
*3 S.K. Ma, C. Dasgupta, and C.-K. Hu, Phys. Rev. Lett. 43, 

1434 (1979). 



